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The energetics of Schottky defects in 123 cuprate superconductor series, REBa2Cu307 (where RE 
= lanthandies) and YAE2CU3O7 (AE — alkali-earths), were found to have unusual relations if one 
considers only the volumetric strain. Our calculations reveal the effect of non-uniform changes of 
interatomic distances within the RE-123 structures, introduced by doping homovalent elements, on 
the Schottky defect formation energy. The energy of formation of Frenkel Pair defects, which is an 
elementary disordering event, in 123 compounds can be substantially altered under both stress and 
chemical doping. Scaling the oxygen-oxygen short-range repulsive parameter using the calculated 
formation energy of Frenkel pair defects, the transition temperature between orthorhombic and 
tetragonal phases is computed by quasi-chemical approximations (QCA). The theoretical results 
illustrate the same trend as the experimental measurements in that the larger the ionic radius 
of RE, the lower the orthorhombic/tetragonal phase transition temperature. This study provides 
strong evidence of the strain effects on order-disorder transition due to oxygens in the CuO chain 
sites. 



I. INTRODUCTION 

It is well known that during the fabrication process 
of superconductor materials, a variety of point defects, 
such as substitution and interstitial impurities, vacan- 
cies, and cation-disorder are involved. These defects have 
a large effect on the properties of superconductors Q . In 
type II superconductors, high critical current densities 
can be achieved by the presence of high-density defects 
which will provide suitable pinning centers for the mag- 
netic flux lines. The ideal size of defects for flux line 
pinning should be comparable to the superconducting 
coherence length. For cuprates such as YBa2Cu307„5 
(Y-123), the coherence lengths are in the order of tens of 
angstroms while the conventional superconductors have 
a coherence length of several thousand angstroms. Thus 
atomic-scale structural inhomogeneities such as point de- 
fects and columnar defects can play an important role in 
flux-line pinning An increasing number of applica- 
tions of the 123-type high- Tc superconductors use mate- 
rials other than Y-123. For example, in many bulk forms 
and multilayer applications, Y is replaced by Nd, Sm, or 
other rare-earth elements (RE). Doping YBa2Cu307_5 
with Ca, Sr, or other alkaline earth elements (AE) has 
also been shown to improve bulk and grain boundary 
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transport and other properties ^ |^ . 

Since it is well-known that strain effects are impor- 
tant in the studies of point defects, it is expected that 
studies of strain effects on point defects for series of 
REBa2Cu307_5 (RE-123) compounds will be important 
for practical applications of superconductivity. In partic- 
ular, the concentration and ordering of oxygen vacancies 
have significant effects on the superconducting proper- 
ties. For Y-123, the superconducting temperature, Tc, 
depends on the oxygen stoichiometry. As an example 
of a generic doping curve in cuprates 0, 0| , when 5 is 
larger than around 0.7, the crystal loses superconductiv- 
ity. However, if 6 is smaller than 0.7, the compound is 
superconducting. It is generally believed that higher oxy- 
gen content can create more holes in the structure. When 
S is between 0.1 and 0.7, the crystal is in the underdoped 
region, and Tc increases with increasing hole concentra- 
tion. When the oxygen content is in the proper range, 
Tc is above the boiling point of liquid nitrogen. However, 
if more holes are created, the Tc value decreases instead, 
and the crystal is in an overdoped region. 

Oxygen ordering in the Cu-0 chains of the 123 struc- 
ture gives rise to further complex structures 0| . Even 
when the average occupancies of oxygen sites remain con- 
stant, the occupation at chain and anti-chain sites can 
vary. Jorgensen et al. ^3 observed that the super- 
conducting transition temperature in Y-123 changes as 
a function of time following the quenching experiment 
while the oxygen content is fixed. This demonstrates 
that the specific ordering of the oxygen atoms in the 
basal plane is another important parameter that controls 
Tc in the Y-123 system ^J. The ordering process in 



2 



the Cu-0 chain is the origin of the structural transition 
between tetragonal and orthorhombic phases in the 123 
structure, which has been extensively studied since the 
discovery of Y-123 [H El El 113 . Structural transition 
of Y-123 is strongly affected by external pressure as well 

[H El El El 

In addition, there are systematic results 
reported by Wong-Ng et al. HElPl showing that the 
orthorhombic/tetragonal phase transition temperatures 
in RE-123 can be scaled approximately linearly with the 
ionic radius of RE'^+. The above experimental observa- 
tions indicate that lattice strain may play an important 
role in phase transition, which has not been investigated 
theoretically and systematically so far. 

In this paper, we plan to study the effects of homo- 
valent substitutions at Y and Ba sites, and hydrostatic 
pressure on the order-disorder transition theoretically. 
Considering large stress fields due to dislocations around 
grain boundaries, this study will also provide valuable in- 
formation for understanding the transport properties in 
the vicinity of grain boundaries. First we briefly explain 
methods used in atomistic simulations. Secondly, we fo- 
cus on the effects of strain on Schottky defects and re- 
lated phenomena for a series of homovalent substitutions 
at Y and Ba sites of YBa2Cu307_5. Although the 123 
structure is not stable in the Ca- and Sr-analogs under 
ambient pressure p3 . l23 . l25l |. partial substitution of Sr on 
the Ba site has been reported 0, i2^ i2^ ,22„ 2^ ■ Finally, 
strain effects on chain-oxygen order-disorder transition of 
REBa2Cu307_5 will be investigated. 



II. ATOMISTIC SIMULATION METHODS 

Since the parent compounds of high- Tc cuprate super- 
conductors can be considered as charge-transfer insula- 
tors, ionic bonding can be assumed to have a large con- 
tribution to the lattice energy. Methods used to study 
atomistic phenonmena in conventional oxides can there- 
fore be applied to study cuprates. The cuprates become 
superconductors at a proper doping level and temper- 
ature. However, the charge carrier density is very low 
compared with that of conventional metals. In addition, 
the charge carriers are confined within copper-oxygen 
planes. Consequently, the screening effect is not as strong 
as that of conventional metals. The lattice energy calcu- 
lated from ionic models is somewhat over-estimated. The 
point defect's energy is also slightly overestimated due to 
the omission of screening effect in metallic region of the 
phase diagram, which can be improved by including po- 
larization effects in the shell model 29^. Many previous 
theoretical investigations are based on this type of ionic 
model (for example,see [H IH HI IH HI ) . 

Several pair-potential sets of shell model parameters 
have been determined for Y-123 by Baetzold [3ll l33| . 
To systematically study homo-valent substitutions on Y- 
and Ba- sites, a consistent set of shell model parame- 
ters using the data set for YBa2Cu307 [^l^was further 
developed to account for the dependence of the Born re- 




FIG. 1: Crystal structure of YBa2Cu307- 



pulsion of the two ions on their net charges, on outer 
electronic configurations 35, 36], and on the common "r'^ 
law" between polarizability and radius |37| . The "virtual 
crystal method" is applied here to interpret experimental 
data of mixing two types of elements on one site. This 
method essentially is, for the purpose of calculating aver- 
age structural and elastic properties, to approximate the 
mixture of two ions distributed over one sublattice by 
identical average "virtual ions". It allows the incorpo- 
ration of compositional changes at the atomic level, but 
ignores explicit effects of disorder. For instance, when 
Ba^^ is partially replaced by Sr^+ in experiments, the 
composition-weighted average value of the two ions' radii 
is taken to approximate that of each virtual divalent ion. 
After constructing a consistent interatomic pair potential 
set 113, we used the "General Utility Lattice Program" 
(GULP) 39], which integrates the above modeling meth- 
ods at an atomistic scale, to study lattice energy, elastic 
constants, and lattice dynamic properties. In summary, 
our calculations are based on short-range potentials of 
the Buckingham type, long-range Coulomb potentials, 
and displacement-induced deformations of the electronic 
charge density in the framework of a shell model. 



III. SCHOTTKY DEFECTS IN 123 
COMPOUNDS 

Any deviation in a crystal from a perfect structure 
is an imperfection. The simplest imperfection is a lat- 
tice vacancy, which is a missing atom or ion, known as 
a Schottky defect. Schottky defects involve "multiple" 
vacancies while preserving electrical neutrality. Regard- 
less of how a Schottky defect is created, it is necessary 
to expend a certain amount of work per atom to take 
it to the surface. We calculated point defect energy by 
Mott-Littleton approach [ifl liH ] . Some vacancies are at 
anionic sites and others at the cationic sites. For Y-213, 
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FIG. 2: Schottky defects formation energy vs volumetric 



strain. The volumetric strain is defined as 



Vo 



'-, where Vo is 



the volume of Y-123 compound. All volumes are obtained by 
optimizing cell parameters and internal coordinates to mini- 
mize the total energy. The solid solutions are treated by the 
"virtual crystal method". 



where c is the equihbrium concentration of vacancies; c is 
the average concentration of vacancies. The equilibrium 
concentration of the vacancy j, Cj, can be computed from 
Boltzmann statistics as 



exp(. 



kT 



(3) 



where Nj is the number of atom j per unit volume, and 
Ey. is the energy required to take atom j from its lat- 
tice site inside the crystal to a site on the surface. Sub- 
stituting the equilibrium concentraions of vacancies into 
equation |(5J), we obtain the expression for the average 
concentration of vacancies, c, as follows: 



c — N ■ exp(- 



-Es chottky 

kf 



); 



(4) 



where N is the number of formula units per unit volume, 
and Eschottky is given, in terms of the formation energy of 
individual vacancy (an ion is removed to infinity instead 
of the surface of crystals) and lattice energy E 
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FIG. 3: Schottky defects formation energy vs Bf2 of 123 
compounds. B is the bulk modulus and is the mean volume 
per atom of 123 compounds. Both B and Q are calculated 
from the optimized structures with minimal total energy. The 
solid solutions are treated by the "virtual crystal method" . 



the defect reaction is given as follows: 



+ 3-V^^ + 7-V^^ ; (1) 



where V represents a vacancy. According to mass action 
law, the equilibrium constant (K) at a finite temperature 
can be written as 



K = c{Vy')-ciV^:^r-c{Vc:f-ciV^^) =c^^ (2) 



13 7 



E, 



Schottky 



'Vacancy 



E 



lattic 



13 



(5) 



In RE-123, the absolute value of site potential de- 
creases continuously with increasing ion radius at Y-site, 
which is consistent with applying external tensile pres- 
sure. Usually, the change of short-range repulsion is less 
significant than that of the Madelung site potential. In 
general, the smaller absolute site potential value is favor- 
able for lowering vacancy energy, and vice versa. Fur- 
thermore, an increase of the total volume leads to the 
expansion of effective relaxation space so that the va- 
cancy energy can further decrease due to the relaxation 
of the atoms. In the shell model, the electronic polar- 
ization due to the electronic relaxation also lowers the 
defect energy. In YAE2CU3O7 (AE-123), as the radius of 
the ion at AE-sites becomes larger, a similar phenomenon 
is observed. In both REBa2Cu307 and YAE2CU3O7, the 
lower Schottky defects formation energies are always as- 
sociated with the larger volume. This indicates that the 
volume of formation of Schottky defects is positive, or 
the volume of the entire crystal expands during the for- 
mation of Schottky defects. 

As shown in Fig. (1), the complex structure of 123 
compounds can be roughly divided into six layers along 
c-axis direction: RE-Cu02(Cu2)-BaO-CuO(Cul)-BaO- 
CUO2-RE. For the sake of simplicity, the distances be- 
tween planes are represented by the separations of cations 
projected along c-axis. The external and internal strains 
are computed with Y-123 as a reference. Replacing Y'^+ 
by larger RE'^+ ions leads to positive external strains, 
indicating the dimensions of the cell expand as increas- 
ing RE size. The more interesting observations are the 
changes of internal parameters within the unit cell. The 
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variance of internal strains provide direct information of 
the changes of layers' separations in RE-series. The inter- 
atomic distances vary quite differently with ionic radius 
in the RE-123 and AE-123 series [s^. For instance, the 
distance of RE-CUO2 becomes larger with increasing RE 
size, which is expected from the change of short-range 
repulsion of RE-0. However, even when the entire cell 
volume increases as RE ion becomes larger, some parts of 
unit cell contract such as interlayer separations between 
Cu02(Cu2)-BaO and BaO-CuO(Cul) (or, bond lengths 
of Cu2-04 and Cul-04). Consequently, the average en- 
ergy of the point defect only decreases slightly because 
of the existing compressive blocks within the structure 
as plotted in Fig. (1). In the AE-123 series, when Ba^+ 
is (partially) replaced by smaller AE^+, the variances of 
external and internal strains are of the same sign, reflect- 
ing a somewhat "even" expansion of the entire structure. 
The difference between RE-123 and AE-123 is reflected 
by the changes of AE-AE and Cul-Cu2 distances. Un- 
like those in RE-123, for AE-123 both distances become 
larger with increasing AE size. The different changes of 
internal strains due to doping at either Y or Ba sites gov- 
erns the absolute value of slope of Schottky defects of RE- 
123 and AE-123 (see Fig. jSJ). In the previous studies 
ji^ . they found interesting correlations between dopant 
radius and energy of solution in YBa^CusOy. The trend 
of subsituting divalent cation for barium matches well 
with that in the Schottky defect energy of YAE2CU3O7. 
However, there exist clear different trends between re- 
placing yttrium by trivalent ions in YBa2Cu307 and that 
in the Schottky defect energy of REBa2Cu307. The rea- 
son is that the parent structure is fixed for former case, so 
that the trivalent ion larger than yttrium leads to more 
energy of solution. One the other hand, the Schottky de- 
fect has to consider all the atoms in the unit cell. Unless 
the structures change more or less "homogeneously" , the 
energy of solution at one specific site may evolve differ- 
ently from the Schottky defect energy. 

Recall that the relation between cell volume, bulk mod- 
ulus of RE-123, and the trend of thermal expansion co- 
efficients are also the manifestation of the complexity of 
the 123 structure [s^. There is an "unusual" relation be- 
tween Schottky defects formation energy and BQ (where 
B is bulk modulus and is the mean volume per atom) 
of RE-123 in Figure ©. As reported by Varotsos [43l |. 
the Schottky defects formation energy is proportional to 
Bft for elemental and binary crystals. While this relation 
appears to be obeyed by the AE-123 series, it is violated 
by the RE-123 series. The possible reason is still the non- 
uniform changes of the interatomic distances within the 
unit cell for RE-123. We have studied the trend of melt- 
ing temperature of RE-123 by "Lindemann law" 0,^3 j 
and found that the vibrations along c-axis are the most 
important modes to determine the melting temperatures 
while the isotropic approximation fails to yield correct re- 
sults ;3Sj. Here, we take another approach to investigate 
the relation between melting temperatures and formation 
energies of Schottky defects. Kurosawa gave a direct 



TABLE L The Schottky defects formation energy and melt- 
ing temperature of binary and RE-123 Compounds. The data 
of binary compounds are from Ref. .46] . The melting tem- 
perature data of RE-123 are from Ref. |47|. 

Crystal NaF NaCl NaBr MgO CaO~SR)" 
Eschottky(eV^) 2.5 2.2 2.0 6.3 5.5 5.0 
Tmciting(°K) 1259 1073 1018 3070 2850 2700 
Crystal YBCO GdBCO EuBCO SmBCO" 
Eschottky(eT/) 4.22 4.21 4.20 4.16 
Tmciting(°K) 1250 1290 1300 1325 



correlation between them, which is listed in table iQJ. 
While it appears that a smaller Schottky defect forma- 
tion energy can lower the melting temperature of simple 
binary significantly, it has not been proved to be true 
for compounds with complex structure. Note that the 
predominant point defect could be Frenkel defect rather 
than Schottky defect in some complex structures. How- 
ever, Schottky defect energy provides an unambiguous 
measure on the average cohesive strength. For example, 
despite the Schottky defect formation energy of Sm-123 
is smaller than that of Y- 123, Sm-123 has a higher melt- 
ing temperature. It is necessary to study internal strains 
rather than volumetric strains in order to obtain a de- 
tailed analysis of the thermodynamic properties of com- 
pounds with complex structures. In general, we found 
that the signs of the slopes (dEschottky/dev) of RE-123 
and AE-123 in Fig. Q are the same. The difference be- 
tween these two slopes and the "unusual" relationships 
between Schottky defects formation energies and Eft and 
melting temperatures of RE-123 compounds all reflect 
the complexity of the crystal structure. 



IV. CHAIN-OXYGEN ORDER-DISORDER 
TRANSITION 

The oxygen content is an important parameter for 
superconductivity in cuprates. In addition, the distri- 
bution (ordering) of oxygen atoms among the atomic 
sites strongly affects Tc- The order-disorder transition of 
chain-oxygen has been extensively studied both experi- 
mentally and theoretically H H El El El 113 . Raman 
studies show that there exists a pressure-induced order- 
ing phenomenon El| ■ The strain effect on this transition 
has not been investigated theoretically in a systematic 
fashion. Here we focus on the effects of homo-valent sub- 
stitutions at Y and Ba sites, and hydrostatic pressure on 
the order-disorder transition. 

In the quasi-chemical approximation (QCA)^l [l5l| . 
the short-range order is characterized by the fraction 
number of near-neighbor pair sites occupied by oxygen- 
oxygen pairs p = , where N is the number of sites on 
each of the sublattice a and /?, and Nqo is the number of 
oxygen-oxygen near neighbor pairs. The long-range order 
parameter S is defined such that the fractional site occu- 
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pancy on sublattice /3 is c{l + S), while that on sublattice 
a is c{l — S) where c is the fractional site occupancy aver- 
aged over both sublattices. [Note that for YBa2Cu307_5, 
c is 0.5 when 5 is zero; i.e. S—l-2c.] 
The partition function is given by 



-W{Ra ,9„„ ) 



(6) 



where Ra — N{1 + S)/4; q^^ is the probability of pairs 
aa; and g{Ra,q^J) is proportional to the total number 
of ways one can divide N entities into four groups of 
aa,a[3^ [3a, and pairs. The configuration energy is 
denoted as W{Ra,q^^)- As usual, we may replace each 
sum by its maximum term in the summation for the sys- 
tem of a great many assemblies. Hence, we have 



(7) 



in which the most probable value q 
by 



of q^^ is determined 



-(lng(i?a,g^ 



W{Ro,,q^J 
kT 



) = 0. 



(8) 



There are three unknowns for a given value of temper- 
ature T and oxygen partial pressure Po2- c, S, p. Three 
equations involving these three unknowns are obtained 
by requiring that: 1) the chemical potential of oxygen 
atoms is the same on both sublattices; 2) the chemical po- 
tential of oxygen atoms is the same in the solid and in the 
gas phase(which consists mostly of diatomic molecules 
but has an equilibrium concentration of atomic oxygen); 
and 3) the free energy of the system is a minimum with 
respect to the fractional number of oxygen-oxygen pairs. 
The followin g de rivations follow the similar procedure as 
described in [43 . Using above conditions, QCA yields 



c{l + S)-p ^ 3 (l + g)[l-c(l-5)] 
^c(l -S)-p^ 4 ^ (1 -S)[l- c(l -I- S)]'' 



ln[( 



l-c(l + 5),3 c{\ + S)-p 



c(l + S) 



-r ( 



p, 



n=ln[(l|^) 



1 



1 

2Ed 



kf - 



1-2C + P ' ' C ' ikT)i 

cil + S)~p){c{l-S)-p). 



kT 



p{l -~2c + p) 
I 



(9) 
(10) 
(11) 



where e is the energy to remove an oxygen atom from 
the gas and place it in the lattice; v is the repulsion 
energy between near-neighbour oxygens; E^; is the dis- 
sociation energy of one oxygen molecule; ^ is equal to 
4.144 X 10^9 Pa(eV)-'^/2^ The desired values of c, S, and 
p are obtained for given values of T and P02 j by simul- 
taneously solving the above equations. We note that the 
order-disorder transition temperature Tod is related to 
the oxygen-oxygen repulsion energy v and the value of 
the average site occupancy c at this temperature by 



kT, 



ln(, 



16c(l 



r); 



(12) 



od •l-4(l-2c)2- 

Frenkel defects involve an atom displaced from its nor- 
mal site into an interstitial site. If the interstitial site is 
chosen as anti-chain site for 123 compounds, this Frenkel 
pair is closely related to chain-oxygen order-disorder 
transition. Forming a Frenkel pair in an otherwise perfect 
crystal is an elementary disordering event. As the disor- 
dering proceeds, it is important to account for defect- 
defect interactions. Using Mott-Littleton approach, we 
have computed the isolated Frenkel pair formation energy 
for REBaaCuaOy, YAE2CU3O7, and Y-123 under hydro- 
static pressure. We found that this formation energy 



increases in compressive regions and decreases in tensile 
regions under hydrostatic pressure. This is the origin of 
ordering under stress. From the systematic investigation 
of the phonon spectral characteristics with the applica- 
tion of pressure it was observed that the changes 
induced by the hydrostatic pressure have a strong effect 
on chain ordering. Results of our calculations are con- 
sistent with the reported observation. Figure Q shows 
the plots of the Frenkel pair formation energy versus lat- 
tice strain (volumetric strain) for both the RE-series and 
the AE-series. These two plots demonstrate a similar 
trend and they form a well-connected smooth curve. The 
similar behavior of the curve for the RE- and AE-series 
(with chemical doping) and the curve for Y-123 (with hy- 
drostatic pressure) indicates that the oxygen disordering 
energy is dominated by lattice strain, which is expected 
if the short range repulsion terms dominate the energy 
required to form the oxygen interstitial ion. 

Furthermore, scaling the oxygen-oxygen short-range 
repulsive energy by the calculated Frenkel Pair formation 
energy, the transition temperature between orthorhom- 
bic and tetragonal phases is computed based on equa- 
tion H12|) . The transition temperature of YBa2Cu307„5 
was used as a reference. The results are plotted in Fig. 
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FIG. 4: Frenkel pair formation energy vs. volumetric strain 
for 123 Compounds. The energy corresponds to moving one 
chain-oxygen into an anti-chain site. All volumes are obtained 
by optimizing cell parameters and internal coordinates to min- 
imize the total energy. The solid solutions are treated by the 
"virtual crystal method". 



FIG. 5: Chain-oxygen order-disorder transition temperature 
of REBa2Cu307. The theoretical data are calculated by using 
the transition temperature of YBa2Cu307_5 as a reference. 
The experimental data are from Refs. |2Ci l2ll . The er- 
ror bars indicate the temperature range measured by X-ray 
Diffraction. 



^ with the experimental data taken from [20l l2ll . 
There exists observable difference (c.a. 100 K) between 
theory and experiment in Nd-123 system It has 

been determined experimentally that the orthorhombic- 
to-tetragonal phase transition in the RE-123 series take 
place at an oxygen composition in the range of 6.4 (Er) 
to 6.83 (Nd), not aU at 6.5 (Y) miUm. The change 
of oxygen content leads to the change of structure such as 
lattice parameters and atomic positions |^ , which in 
turn alters the formation energy of Frenkel pairs. Since 
we set the reference transition temperature to be the 
value in Y-123 case, this corresponds to the transition 
at oxygen content being 6.5. If we track the whole pro- 
cess of order-disorder transition starting from fully oxy- 
genated case, the total oxygen content decreases until the 
transition is finished. The formation energy of Frenkel 
pairs decreases also as anisotropy in the ab plane (de- 
fined as (b — a)/a) reduces. Note that the transition 
occurs at oxygen content being 6.83 in Nd-123 system, 
this indicates the anisotropy in Nd-123 remains higher 
compared with Y-123. Therefore, the formation energy 
of Frenkel pairs decreases less than that in Y-123 during 
the whole process of transition as observed experimently. 
This means that in our "simple" model the transition 
temperature for Nd-123 is underestimated. But, the re- 
markable thing is that using such a simple model, the 
trend of the theoretical results agree well with the exper- 
imental measurements. It is seen that stress and "chem- 
ical pressure" can substantially alter the degree of disor- 
der. This study provides a clear evidence for the effects of 
strain on order-disorder transition. Previous studies also 
show there are quite rich microstructures resulting from 
this orthorhombic-to-tetragonal phase transition: for in- 
stance twin structures and related twinning dislocations 



^^113, tweed morphologycaused mainly by (110) and 
(110) shear displacements [l|. In particularly, our results 
can be applied further to deduce that strain near edge 
dislocations in low angle grain boundaries [sij . which 
can also have a significant effect on the degree of order- 
ing in cuprate materials. For example, the interactions 
between point defects and strain fields due to disloca- 
tions and/or grain boundaries can affect the distribution 
of point defects, and the content and degree of order of 
oxygen atoms [3^. 

V. CONCLUSIONS 

Based on the Mott-Littleton approach we studied the 
Schottky defect formation energy in the 123 phase as 
a function of volumetric strain. Generally, a more ex- 
panded lattice favors a lower Schottky defect forma- 
tion energy, and vice verse. The difference of slopes 
(dEschottky/de) between RE-123 and YAE2Cu307_5, the 
"unusual" relation between Schottky defects formation 
energies and and melt temperatures of RE-123 com- 
pounds all reflect the complexity of the crystal structure 
of 123 compounds. 

Our study also illustrates the importance of strain ef- 
fects on the orthorhombic/tetragonal phase transition in 
the RE-123 compounds. We have calculated the forma- 
tion energy of Frenkel pair defects as a function of vol- 
umetric strain for REBa2Cu307 and YAE2CU3O7, and 
for YBa2Cu307 under hydrostatic pressure. Our calcu- 
lations show good agreement with experimental observa- 
tions in that pressure favors ordering of the CuO chains. 
For example, the Frenkel pair formation energy indeed 
increases significantly (c.a. -0.25 eV / 0.01 volumetric 
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strain) under compression. Based on a quasi-chemical 
approach, the orthorhombic/tetragonal transition tem- 
peratures for RE- 123 have been computed by scahng the 
effective oxygen-oxygen short-range repulsive energy in 
the CuO chain using the Frenkel pair formation energy. 
The calculated results agree with experimental data in 
that the larger the ionic size of RE, the lower the or- 
thorhombic/tetragonal phase transition temperature. 
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